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Abstract. Using numerical simulations we show that a variety of novel colloidal 
crystalline states and multi-step melting phenomena occur on square and 
triangular two-dimensional periodic substrates. At half-integer fillings different 
kinds of frustration effects can be realized. A two-step melting transition can 
occur in which individual colloidal molecules initially rotate, destroying the overall 
orientational order, followed by the onset of interwell colloidal hopping, in good 
agreement with recent experiments. 
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Colloidal particles are an ideal system for studying 2D ordering and melting for 
different kinds of substrates as the individual particle positions and dynamics can 
be directly visualized. Crystallization and melting on 2D periodic substrates is also 
relevant to vortices in superconductors with periodic pinning arrays [0, ||, || and atomic 
orderings jjj. Colloidal crystallization on 2D periodic substrates has been the subject 
of considerable recent interest. In several recent experimental studies, a 2D substrate 
for colloids was created using optical tweezer arrays jH), templating || |t|, and 2D 
crossed laser arrays ||. Colloidal ordering and melting on 2D periodic square and 
triangular substrates has been studied through simulation Q and experiment ]Io| . 
It was shown that a rich variety of novel colloidal crystalline states, referred to as 
colloidal molecular crystals (CMC's), form at integer filling of the periodic substrates. 

In Ref . H , integer matching states up through a filling of four colloids per minima 
were illustrated for both square and triangular arrays, and the two-stage melting 
behavior of the dimer state at a filling of two colloids per minima was demonstrated 
for the square substrate. Here, we present the structures observed at half-integer 
filling for the square and triangular substrates, showing that the perfect regularity of 
the integer filling state is destroyed away from integer filling. We also demonstrate 
a two-stage melting transition for the trimer state at a filling of three colloids per 
minima for the square substrate. 

We simulate a 2D system of N c colloids with periodic boundary conditions in 
the x and y directions, using Langevin dynamics as employed in previous colloidal 
simulations fill . The overdamped equation of motion for a colloid i is 
(hi 
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Here = — Y^iJj ^iV( r ij ) is the interaction force from the other colloids. 
The colloid-colloid interaction is a Yukawa or screened Coulomb potential, V(rij) = 
{Q 2 /\?i — i"j |) exp(— «;|rj — rj\), where Q = 1 is the charge of the particles, 1/k is 
the screening length, and ^Uj) is the position of particle i (J). The system length 
is measured in units of the lattice constant ao and we take the screening length 
1/k — ao/2. For the force from the 2D substrate, we consider both square and 
triangular substrates with strength A, period ao, and N m minima. For square 
substrates, f s — Asm(2Trx/ao)x + Asin(27ry/ao)y, and for triangular substrates, 
f s = J2i=i 3Asin(27rpi/ao)[cos(0j)x-sin(0i)y], where = xcos(8i) -ysin(^) + a /2, 
9i = 7r/6, 62 = 7r/2, and #3 = 57r/6. The thermal force fr is a randomly fluctuating 
force from random kicks. We start the system at a temperature where all the colloids 
arc diffusing rapidly and gradually cool to T = 0.0. We do not take into account 
hydrodynamic effects or possible long-range attractions between colloids. 

The colloidal positions for a system with a square substrate at integer matching 
between the colloidal periodicity and the substrate periodicity, N c = nN m , are 
illustrated in Fig. 1 of Ref. || for n — 1 to 4. In particular, at N c = N m , each 
colloid is located at the center of the potential minima and a square colloidal crystal 
forms, while at N c — 2N m , each minima captures two colloids which can be regarded 
as a colloidal dimer with a rotational degree of freedom. Over a range of substrate 
strengths, the colloidal dimers form a rotationally ordered state, with neighboring 
dimers perpendicular to one another. The orientational ordering of the dimers is due 
to the colloidal repulsion, and allows the distance between the colloids to be maximized 
under the constraint of the square substrate. In Fig. 1, we show the colloidal positions 
for the square substrate system at a half-integer filling, N c — l,5N m . Here we find a 
combination of the single colloids and dimer colloids. The single colloids are located 
at every other site with intervening dimers, creating a checkerboard ordering. The 
strong orientational ordering of the dimers observed at N c = 2N m is lost in this state. 
The dimers alternate their orientation between vertical and horizontal in a disordered 
fashion, forming a pattern of grain boundaries. 

The ordered colloidal crystalline states that form on a triangular substrate at 
integer matching from n = 1 to 4 are illustrated in Fig. 2 of Ref. || . As in the case 
of the square substrate, at N c = N m each minima captures a single colloid, while at 
N c = 2N m each minima captures two colloids which form a dimer state. The dimers 
again have an additional orientational ordering in which the dimers in each row have 
the same orientation, which is rotated 45° with respect to the adjacent rows. In Fig. 
2 we illustrate the colloidal positions for the triangular substrate at half-integer filling, 
N c — 1.5N m . Unlike the square substrate, where dimers could tile the lattice in an 
orderly way, filling every other minima, for the triangular substrate the positioning of 
the dimers is frustrated. It is not possible to arrange a state such that minima with 
only one colloid will be between every other dimer. As a result, the colloids form a 
very disordered state in which the dimers show a wide range of orientations, and not 
merely two as in the case of the square substrate at N c = 1.5N m , or the triangular 
substrate at N c — 2N m . 

In Fig. 3 we show the two-stage melting of the CMC at N c = 3N m on the square 
lattice. As illustrated in Fig. 3(a) at low temperatures T/T^ < 0.25 (where is the 
melting temperature at zero substrate strength), both orientational and translational 
order of the trimers are present and the system is frozen. This is the "ordered solid" 
phase. In Fig. 3(b), at T/T^ = 1.5, the trimers begin to rotate within the minima; 
however, diffusion of individual colloids throughout the sample does not occur. The 
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system is still frozen but the trimer orientational order is lost. This is the "partially 
ordered solid" phase. In Fig. 3(c), for a higher temperature T/T^, the system enters a 
modulated liquid phase. Here the colloids begin to diffuse throughout the system. We 
note that in the recent experiments of Brunner and Bechinger fllPf , the same melting 
phenomena was also found for trimers on a triangular substrate. 

To summarize, we have shown the rich variety of novel colloidal crystalline and 
disordered states that can be achieved with square and triangular two-dimensional 
substrates. The colloidal molecular crystal states that appear at integer filling of the 
substrate minima persist at half-integer filling, appearing as a mixture of the two 
neighboring integer filling states. For the case of N c = 1.5N m , the single and dimer 
states coexist. An ordered filling of dimers can be arranged on the square substrate at 
N c = l.5N m , but the dimer arrangement is frustrated on the triangular substrate and 
a strongly disordered state results. The triangular lattice at N c = 1.5N m is thus a 
realization of a frustrated system which would provide an interesting topic for further 
study. For integer filling of N c = 3N m on a square substrate, we demonstrate the 
multistage melting of the CMC, where the orientational order of the colloidal molecule 
states is lost first, followed by the translational order. Since the colloids within a 
minima can act as a single particle with a rotational degree of freedom, our results also 
suggest that certain canonical statistical mechanics models, such as Ising, XY, Potts, 
and frustrated models, may be realized with colloids on two-dimensional periodic 
substrates. The states predicted here should be observable for colloids interacting 
with crossed-laser arrays or optical tweezer arrays, dusty plasmas in 2D with periodic 
potentials, and vortices in superconductors with periodic substrates. 
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Figure 1. The colloid configurations (black dots) at T = 0.0 for a square 2D 
periodic substrate with A = 2.5, for a half-integer colloidal density, N c = 1.5N m - 
Every other minima captures a single colloid, while the remaining minima capture 
two colloids in a dimcr state. There is no long-range rotational order of the dimcrs. 
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Figure 2. The colloid configurations (black dots) at T = 0.0 for a triangular 2D 
periodic substrate with A = 2.5, for a half-integer colloidal density, N c = 1.5N m . 
It is not possible to arrange a state such that minima with only one colloid will 
be between every other dimer. Therefore there is no dimcr ordering. 
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Figure 3. The colloid positions (black dots) and trajectories (lines) over fixed 
time intervals at different temneratures for the trimer state on a square substrate 
shown in Fig. 1(c) of Rcf. ||. (a) The ordered solid phase at T/T^ = 0.25, 
where is the melting temperature at zero substrate strength, (b) T/T^ = 1.5. 
Oricntational order is destroyed as the trimers rotate within the substrate minima, 
but the colloids remain trapped inside each minima so the system is still in a 
solid phase, (c) T/T^ = 4.0. Individual colloidal diffusion occurs throughout the 
sample and the system is in the liquid phase. 



